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I. Introduction 



In this paper we consider the nonrelativistic many-electron system in the forward, 
exchange and BCS approximation. In this approximation, which is still quartic in the an- 
nihilation and creation operators, the model can be solved explicitly. The partition func- 
tion and the correlation functions are given by finite-dimensional integral representations. 
We work in the quantum grand canonical ensemble and start with positiv temperature 
T = jj > and finite volume L d < oo. 

The standard model of the many-electron system in d space dimensions is given by 
the Hamiltonian 

H = H + H int (1.1) 

where 

#o = -h Yl ek a L a k* (1-2) 

k,cr 

and 

H ijlt = jjz Yl ^( k -P) a k«7 a P* a J-kT a q-PT ( L3 ) 

o-,t€{T,I} k,p,q 

We assume U to be short range, that is U is L 1 in coordinate space. The energy momentum 
relation ek is given by 

e k = |^-M (1.4) 

or may be substituted by a more general expression which satisfies ek = e_k- The param- 
eter ji > 0, the chemical potential, is present since we are working in the grand canonical 
ensemble and is determined by the density of the system. Since we are in finite volume, the 
spatial momenta k range over some subset of (^Z) . The physics of the nonrelativistic 
many-electron system is determined by momenta close to the Fermi surface 

F={ke (fZ) d |e k = 0} (1.5) 

so we impose a fixed ultraviolet cuttoff and choose 

kGM w = {k£ (ir Z ) d | l e kl (1-6) 

In the context of conventional superconductivity the cuttoff ui is referred to as the Debye 
frequency. 

The normalizations choosen in the definition of the creation and annihilation operators 
are such that the anticommutation relations read 

{a k <7, a£> a <} = £ d £k,k'£<r,<r' (1-7) 
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In particular, if Op = Ua^ a 1 is the zero temperature ground state of the noninteracting 

krr 
e k <0 

system, then (Op, a^ako-Op) = L d 9{— ek) such that (Op,i?oOp) = ^ k ct ek 0(— ek) ~ 
ek ^(— ek) is indeed an extensive quantity as it should be. Here 6(e) is the step 
function being one for e > and zero otherwise. 

We are interested in the grand canonical partition function 



Z = Z(P,L,U) = Tre~ pH 



(1.8) 



which may be normalized by Tr e such that Z(U = 0) = 1, and in particular in the 
two point function 

Tre-Wa^avv, 



(a k( r a k'(T')/?,L 



(1.9) 



Tre-I 311 

which gives the momentum distribution of the system. Recall that in the free (U = 0) 
system, the ideal Fermi gas, 



lim (a£ ' ak' a ')p jL = L Su,k' S a,a'0(-eu) 



(1.10) 



Since we are in the quantum grand canonical ensemble, the traces in (1.8,9) are to be 
taken over the Fock space T = ©^L -^n where 



•F'. n 



|F n e L 2 ([[0,L] d x {T,|}] n ) Fn^KiTr!, ■ ■ ■ ,-XL nn a nn ) = sgn7rF n (xicri,---,x n cr n )| 



The fact that the physical system has a fixed number of particles iV is expressed by 
requireing that the expectation value (N)^ p of the number operator N = J^kcr a k<T°kc r 
(which is extensive, see above) is equal to the number of particles, (N)^ = N, which 



determines \x as a function of the density 

As usual, the quartic part H int (1.3) of the Hamiltonian H (1.1) may be represented 
by the following four legged diagram: 



k, a 



q-k,T 



int 



C/(k-p) 



P,(T 



q- P,r 



Because of conservation of momentum, there are three independent momenta here labelled 
with k, p and q. Then one can consider the following three limiting cases with only two 
independ momenta: 
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forward 



exchange 



BCS 



k, a 



k, a 



k, a 



-k,T 



U(0) 



C/(k-p) 



t/(k-p) 



k, a 



P,(T 



k,r 



P,(T 



-P,T 



That is, one may consider the approximation 



Hint ~ -ffforw + ^ex + -f^BCS 



Hbcs 



(1.11) 



where 



#forw = 5^ E ^ ~ P ^ k 

ct,t6{T,4} k,P,q 

= 1^ X E a k>kaa+ r a pT 



P °kcr a P^ a q-kT a q-pT 



Y ~ P) 5k ><l-P a ka a P- a q"-kr a q- 

ct,t6{T,4} k,p,q 

1] J] ^ ~ P ) a k<rV a Jr«kr 



pr 



°Yr€{T,|} k,p 



#BCS = JJI Y ~~ P ) a ka a P^ a q"-kr a q- 

o-,t€{T,I} k,p,q 

o-,t€{T,J.} k,p 



pr 



a kcr a P^ a -kr a -p-r 



(1.12) 



(1.13) 



(1.14) 



Let us shortly make a comment on the volume factors. In (1.12-14), we only introduced 
some Kroenecker delta's but we did not cancel a volume factor L d . That this is the right 
thing to do, that is, that the left hand side as well as the right hand side of (1. 11) is indeed 
proportional to the volume may be seen in the easiest way for the forward term. On a 
fixed n particle space T n the interacting part H int is a multiplication operator given by 



-ftintl 



= \ Y u (*i-xj) 



(1.15) 



i,3 = l 



Let <5 y (x) = 8(x — y). Then <^(xi, • • • , x n ) = S yi A • • • A<5 yii (xi, • • • , x n ) is an eigenfunction 
of ifi n t with eigenvalue 

n n 

s = SE ^ - yi) = *h E E el(yi " yj)q ^ ( L16 ) 

— l <,j'=l q 

which is, for U E L 1 , proportional to n or to the volume L d for constant density. One 
finds that <p is also an eigenvector of the forward term, Hf orw Lp = E{ orw Lp where Ef orw is 
obtained from (1.16) by putting q = without cancelling a volume factor, 

n 

£forw = ^ E t/ ( C l= ) ( L17 ) 

i,J = l 

which is also proportional to the volume. 

We now come to the exact definition of the model which is solved in this paper. To 
do so, we need the functional integral representation of the perturbation series for the 
partition function. It is summarized in the following theorem which is fairly standard. 
One may look in [FKT1] for a nice and clean proof. 

Theorem 1.1: Let H = H + H int be the Hamiltonian (1.1-3), let 

Z = Z(/3, L) = Tr e-W+Hint) / Tr e -(3H ^ ^ 

be the normalized grand canonical partition function and let 

a) The partition function Z has the following perturbation series 

Z = J2 nV E E II { S k 2i - 1+ k 2i , P2i - 1+ p 2i U{k 2i - P2i )} X 

71=0 <J\"-<J 2 n fe l" ,fc 2ra i=l 

detiS^S^Cih)}^..^ (1.20) 
Here k = (k , k) e j|(2Z + 1) x { k e (^Z) d | \e^\ < u } and the covariance C is given by 

C(k ,k) = — (1.21) 

zfc - e k 

// ||^(x)|| L i < c{f3L d )~ 1 or ||t/(x)|| L oc < c([3L d )- 2 , c a constant, then (1.20) converges. 
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b) Let 

ducty, = n e-& E k J^-e^ k ^ n ^ ^ (L22) 

k,a k,a 

be the Grassmann Gaussian measure with covariance C . Then the perturbation series 
(1.20) can be rewritten as 

Z = J e-v^W^dncty, (1-23) 

where 

°ki+k 2 ,p 1 +p 2 U(k 2 ~ P2)^k 1 a 1 ^k 2 a 2 ^p 1 a 1 ^p 2 a 2 (1-24) 

0\0 2 k 1 ,k 2 
PI .P2 

c) The momentum distribution (1.19) at temperature T = ^ is given by 

{at a a\L*)l3,L = lim I ^ {^k ^k ^a) p,L,e (1-25) 

<=<o fc 6f(2Z+l) 

w/iere 

Here C e (k) = and Z e = J e ~^d^. 

Remarks: (i) A bound from which convergence of the perturbation series follows for 
sufficiently small U is easiest obtained in coordinate space. There the perturbation series 
reads 

z = y e -/«ide a (^)«i)f«i-€ a )(^)(&) d/Xc ( V ,^) 

oo „ n 

= EA / ^•■■^II C/ fe.-i-6i) det[C(ei,e,)] 1 < iJ < 2n (1-27) 

n=0 i=l 

where £ = (x ,x,a), f d£ = J2ae{hl} I[o,/3] dx o I[o,L] d ^ U ^ = U{x) = 8(x )U(x), 
C{£,Z') = &*,<t'C{x-x?) where 

C7(x , x) = Jj ^ e*^-**" [0(-x o )0/j(-e k ) - (9(x )^(e k )] (1.28) 

and 6*^ is an approximate step function given by ^(e) = 1+e 1 _ )gE . First observe that, 
since sup Xo6[0)/3 ] |e" ekX ° [^(-x )^(-e k ) - 0(x o )e f3 (e k )} \ < 1 for all k, one has ||C||oo < 
XlkeM^ 1 < oo if we choose a fixed ultraviolet cuttoff in (1.6). The determinant can be 



bounded using | det[oi • • -a r ] \ < Hl =1 1 1 ^ 1 1 2 where the di are r component vectors. That 
is 



det[C(&,0)] 2n < if { £ \C(^Q\ 2 Y < n(2n||C||L)' = (2n||C||L) n (1-29) 
which gives 

00 00 
\Z\ < E (^) n (4/3^) n |l^ll2 1 (2n||C|| 2 0O ) n = EN^II^T W U W^ (1-30) 



n=0 n=0 

(ii) The e-limit shows up in (1.25) because of 

C(x ,x) IO = /0 ^^e !k V« [^(-e k )(9(-x ) - ^(e k )(9(x )] 

k 

= e *x-*fco*o i_ ( L31 ) 

k,fc 

but 

C(0,x) = £ ^e l ^(-e k ) = lim _^^ e " ()e _^_ (L32) 

k «<o k,fc 

see the proof of Lemma Al in the appendix for a more detailed discussion. 

(iii) Whereas the momenta in the Hamiltonian are d dimensional, the variables in the 
perturbation series are d + 1 dimensional. The additional ko variables in (1.20) are the 
Fourier transform of the x$ variables in (1.27) which in turn enter the perturbation series 
because of 

^ e -0(ffo+Affi„ t ) = ^ e - x ° H °H int e-^~ x ^dx , (1.33) 
that is, because H[ nt does not commute with Hq. 

As remarked under (iii), the momenta in 

#int = X^ E E U ( k -P) a L a P<r a i-kT a <i-PT (1.3) 
cr,r6{T,4} k,p,q 

are d dimensional, but the variables in 

Vint = (j9 ^d)3 E E ~~ P)^ka^pa^ q -kT^ q -pr (1-34) 



are d + 1 dimensional. On the Hamiltonian level, the forward, exchange and BCS approx- 
imation (1. 11) does not lead to an explicitly solvable model, but on the functional integral 
level it does. This is the main result of this paper. 



Theorem 1.2: Let 

{Vf 0r w+V ex + VBCS}(^, = ^2^2 U ( k ~P) [ S k,p + $k,q-p + 8q,o] $ ] pA ] q-kA 'q-pr 

or kpq 

(1.35) 

Then the approximation 

V int « V forw + V ex + VBCS (1-36) 

in (1.23,26) makes the model explicitly solvable. The partition function, the two point func- 
tions and the generating functional for the connected amputated Greens functions are given 
by finite- dimensional integral representations. The dimension of these representations is 
proportional to J if the electron electron interaction is given by 



U(k-p) = { 



e=-j 
j i 



=■. Y,^yi( k )yi(p) (i-37) 



J2 E A,F £m (k')y £m (p') zfd = 3 



1=0 

A£ i em y*. ) i irn yp ) ij u — o 

=0 m=-i 



Here k = |k|(cos<£k,sin<^k) for d = 2 and k' = k/|k| for d = 3. 

The general solution is written down in Theorem II. 1. In the following, we summarize 
the results for Vint ~ Vbcs and for an electron-electron interaction given by (1.37). This 
case is treated in section III.l. 

In that case the two point function (1.26) becomes (we suppress the e): 



/ Fpi^e-^Wu II du l aT dv l aT 
(Mp>l)p,L = PL%, P f^=° (1.38) 

J e -f3L*VM U n du l aT dv l aT 



<JT 1 = 



Here <p l aT = u l ar + iv l ar , ar G {TT, U, Tl}, 



— ( i Po+ e p)[Po+ e p+ <I> -PTx^-pn+^Pii $ Pix] / t qq\ 

P W {Pl+el+n+)(pl+el+n-) V"™> 



J J 

i , _ < — > i 



*k U = E A ^W k ')' ^u=E A ^u^( k ')' ( L4 °) 



1=0 1=0 



<IW = Y,*UL \VW) - Mi-*)} > **** = E [^(k') - 2/K-k')] (1.41) 



2=0 



1=0 



and are the solutions of the quadratic equation 

O 2 - ($ kU $ kU + $_ kU $_ kn + $ kTT $ kTT + $ kU $ kU )o + $ kTT $ kU $ kU $_ kn 

+ $ kTT $ kU $ kU $- kn + $ kTT $ kU $ kTT $ kU + $ kU $ kU $_ kU $_ kU = (1.42) 

Observe that $ kcrT is not necessarily the complex conjugate of $ k(T r, depending on the 
signs of the coupling constants A;. The effective potential Vp is given by 



w) = X>'„i 2 + Kti 2 + K,i 2 )- £ / ws 1 ^ 



cosh(|^g+^) 
cosh ^ej,. 



(1.43) 



For a pure even interaction, that is, if A^ = for all odd angular momentum £, one 
has, if yi(— k) = (— 1) yj(k), 3»ktT = = and the expectation value simplifies to 



/7 / \ VTd x - Po+ e l+*PU*pU 



z=o 



f e-PL'Ve&n) U du l n dv l n 



(1.44) 



with an effective potential 



cosh( § ^/ e^+$ kn $ k ~ ) 



cosh 



(1.45) 



Recall that the free two point function is given by (V^o-Vy ,°) P, L = 0L d 5 p yl/(ipo — e p ) = 
—(3L d 8 PyP i{ipQ + e p )/(pQ + ep). In particular, for a delta function interaction one obtains 
the two dimensional integral representation 



(^p,a^p',a)/3,L = -f3L d 5 P;P ' 



J Po+ e i>+ x ( u + v ) 



| e -/3Ld Vl3 (u,v) dudv 



(1.46) 



where 



log 



cosh(| v /e^+A(u 2 +^ 2 )) 
cosh |ek 



(1.47) 



A positive A corresponds to an attractive interaction. 
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Since the only place where the volume L shows up in the integral representations 
(1.44,46) is the prefactor in the exponential e~@ L the computation of the infinite 

volume limit comes down to the determination of the global minimum of the real part of 
the effective potential V as a function of the (f> l aT : s. 



figure 1 



For repulsive coupling A < Vp (1.47) has a global minimum at u = v = which 
results in 

e -l3L d V {u,v)j J 

lim r -— = S(u)S(v) (1.48) 

and 

lim ^v^V,a)^L = -pS Po y o S(p - p') f£g? (1.49) 

For attractive coupling A > and sufficiently small T = jj the effective potential (1.47) 
has the form of a mexican hat. This results in 

lim (Vw</VV>/3,l = -P6 p0iP , o 5(p - p') p i P + ° e tl\, (1-50) 



p 



where A 2 = Ap 2 , ~ e * and Pq = Uq + Vq is the value where Vp takes its global minimum 
which lies on a circle in the u, v plane. 

The expectation values (ipip) and (ipip) can also be computed. To make them nonzero, 
we introduce a small external field r = |r|e m . That is, we substitute Vbcs by (we do not 
consider here (a^a-^) expectations) 

V B cs,r = Vbcs + jfj ^[rip^ip-ki - fipk^-ki] (1-51) 

k 

One obtains again a finite-dimensional integral representation. For a delta function inter- 
action, that is for J = in (1.37), it reduces to a two dimensional integral: 

/ e-P Ldy e>^dudv 

(M-P',1)P,L = PL%y P ° £P J e U -f3L«V, A u,v) dudv (1-52) 
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where 



cosh( § y/el+\(u 2 +v 2 )) 
cosh 



(1.53) 



Consider an attractive A > and sufficiently small T = jj such that the global 
minimum of the effective potential moves away from zero. For r = 0, the global minimum 
of Vp jr is degenerated and lies on a circle in the u, v plane. In particular, Vg o is an even 
function of u and v and (?/y|-^_ p j) vanishes by symmetry. 



figure 2 



For r 7^ 0, Vg jr has a unique global minimum at (u, v) = (0, t>o) where vq is given by 
the negative solution of 

which, in the limit r — > 0, becomes the BCS equation for |A| 2 = Xvq. Thus 

hm hm - — —j- — - — = hm lim - Hjd ,., — : — r— = — — ri = o(w)d t> + ^ 

(1-55) 

and (ipp^ip-pi) becomes nonzero. 

For repulsive A < 0, Vg >r is complex and the real part of Vp^ r has a unique global 

minimum at (u, v) = (0, 0) which results in lim| r |^o lim^^oo -z — zraZ^v — 17TV) — 

— = S(u)S(v) 

and lim| r |^olimL^oo(^pT ? /'-p4) = °- 

These results of course are also obtained if one applies the usual mean field formalism 
which is based on approximating the quartic Hamiltonian by a quadratic expression. How- 
ever, the situation is different if the electron-electron interaction contains higher angular 
momentum terms and the space dimension is 3. In that case the standard mean field 
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formalism [AB,BW] predicts an angle dependent gap but, if e k has 50(3) symmetry, this 
is not the case. 

Suppose that > is attractive and 

i 

U(k -p) = \ e Y em (k')Y em (p') (1.56) 

m=-£ 

Then the Anderson Balian Werthammer mean field formalism gives 

^lim Ld (a k(T a kfT ) - 2 ^1 e k ^ V<+KAk \ aa ) 

where the 2x2 matrix A k , A k = — A_ k , is a solution of the gap equation 

A p =/ ^U( P ' - k QA k tanh( ^ e ^ Ak) (1.58) 

In 3 dimensions, it has been proven [FKT2] that for all I > 2 (IV. 58) does not have unitary 
isotropic (A k A k = const Id) solutions. That is, the gap in (1.57) is angle dependent. 
However in Theorem III. 3 it is shown that, for even i, 

L lim ^{M^l = kl+<+x ^^ LY , m{xW ^ (1-59) 

which gives, using (1.25), 

if A(x) = A|po S m CKm^m(x) and po and «° are values at the global minimum. The point 
is that for 50(3) symmetric e k also the effective potential has £0(3) symmetry which 
means that also the global minimum has £0(3) symmetry. Since in the infinite volume 
limit the integration variables are forced to be at the global minimum, the integral over 
the sphere in (1.59,60) is the averaging over all global minima. 

There is some physics literature [B,BZT,BR,G,Ha,W] which investigates the relation 
between the reduced but still quartic and not solvable BCS Hamiltonian 

#BCS = #0 + -jjd Yl Yl ~ P ) a k, a -kr a P^-pr ( L61 ) 

°-,Te{i\i} k, P 

and the quadratic, explicitly diagonalizable mean field Hamiltonian 

H M f = H + j^ Yl $^^( k - P) (a+* a -kT( a p° a -pr) + (a kfJ ai kr )a Pf Ta- pr 
°yr€{T,4} k,p 

-( a k, a -kr)(V a -pT}) ( L62tt ) 
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where the numbers (a pcr a_ pr ) have to be determined by the condition 



{q>p<jCI— p T ) — 



jy e -f3H MF 



(1.626) 



which is equivalent to the gap equation (1.58) if one defines A p(TT = -jj ^ (P — k) -jj ( a kcrO- 
One has, if H' := Hbcs — Hmf, 



Opa^—pr 



(a p<T a_ pr )j (1.63) 



It is claimed that, in the infinite volume limit, the correlation functions of both models 
should coincide. More precisely, it is claimed that 

lim log y r e e "-Z B M c F £ 



(1.64) 



vanishes. To this end it is argued that each order of perturbation theory (with respect 
to H') of Tr e -^( H MF+H')/7y e -/3H MF i s finite as the volume goes to infinity. The Haag 
paper argues that spatial averages of field operators like _/j L ^ d <i d xa|(x)aj(x) may be 
substituted by numbers in the infinite volume limit, since commutators with them have an 
extra one over volume factor, but there is no rigorous controll of the error. At least for a 
more complicated electron-electron interaction (1.56), this reasoning cannot be correct in 
view of (1.57-60). 

Namely, consider the (a^aka) p,L expectation. In terms of Grassmann variables it is 
given by (1.25,26). Assume first a delta function interaction U(k — p) = A. For the full 
model (1.3), by making a Hubbard Stratonovich transformation, (1.26) can be rewritten as 



where the integrand is given by 



-V( 



n d(f) qoCl d(f) qoq 
90 q 



fe- v M n# 9oq # 



(1.65) 



90 q 



90 q 



F(d>) = 



and the effective potential reads 



d 


det 


A + S l 


det 


}ijt*) h 4> A 





(1.66) 



V(<t>) 



E i^'/^q 
go q 



109 



-log 



det 


}{^)^4> A 


det 


'A 0" 




A 





(1.67) 
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Here A and S a are diagonal matrices with entries = iko — e k and Sk<? respectively and 
is a short notation for the matrix {4>k-p)k,p where k = (k , k) e |j(2Z + 1) x M u . Of 

course, for example det ^ = Hk(ko + e k ) does not make sense, but the quotients in 
(1.66,67) are well defined. 

The correlation functions of the reduced, but still quartic BCS Hamiltonian (1.61) 
are obtained from (1.65-67) by first assuming that the global minimum 0™ in of V(4>) is 
proportional to 8q,,o4>qo an d second by suppressing the quantum fluctuations around — 
for all q ^ 0. That is, for the model (1.61) one finds 



fF p (j>) 



-L d VW 



Ud(p qo d(p qo 
qo 



Je- Ldv Wnd<f> qo d$ 



Qo 



Qo 



(1.68) 



where the integrand is given by 



F p (<t>) 



and the effective potential reads 



fl 9 det 

O8 vo" |s=0 


~A P + S pr " 


det 


A P i(^<t>' 

}{^) h 4> A P 





(1.69) 



n0) = Ei^oi 2 -iE lo s 



90 



det 






det 





" 

Ak_ 





(1.70) 



The volume factor L~2 in the determinant in (1.66,67) has been transformed away by a 
substitution of variables such that it shows up in the exponent in front of the effective 
potential in (1.68). The matrices in (1.69,70) are labelled only by the ko,po variables, that 
is, A k is the diagonal matrix with entries A^t~ = iko — and is a short notation for 
the matrix {4>k -p )k ,p - 

Contrary to the full model, the volume dependence of the model (1.61) or (1.68-70) is 
such that in the infinite volume limit the integration variables 4> qo in (1.68) are forced to 
be at the global minimum of (1.70). 

The model discussed in this paper (that is, only the BCS part) is obtained from (1.68- 
70) by assuming that the global minimum of (1.70) is proportional to S q0}0 . In that case, 
the only integration variable which is left in (1.68) is the q = mode <j> = 4> and the 
expressions (1.68-70) reduce to the integral representation (1.46,47). 
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Now assume that the elctron electron interaction is given by (1.56) and suppose for 
simplicity that £ is even which suppresses tttt an d ill! contributions in (1.3) and (1.61). 
In that case the model (1.61) gives 



-L d V{<j>) 



■pa I 



where the integrand is given by 



FM = 



and the effective potential reads 



n n#™#™ 

m= — t qo 



j e -L*v{4>) n n#;# 



m=— I qo 



m 
Qo 



d 

9s po° \s=0 


det 


~A P + S pr i($)*$ p " 
_i(|)** P A p + S pl _ 


det 







(1.71) 



(1.72) 



m=—i qo k 



det 


A k i(|)*$k" 
i(^)*$k A k 


det 


"A k " 
A k _ 





(1.73) 



Here $ k denotes the matrix with entries $k, Po po = £m=-^™ -j,'^Tn(k), labelled by the 
Po,p'o variables. The model discussed in this paper is obtained from (1.71-73) by assuming 
that the global minimum of (1.73) has only nonzero </>^ =0 modes. In that case (1.71-73) 
reduce to the integral representation (1.44,45) which further can be reduced, by Theorem 
III. 3, to (1.59). However, the argument used in Theorem III. 3 that the (^pcrV'po-) expecta- 
tion has to be SO(3) invariant if e k is SO(3) invariant still applies to (1.71), since (1.73) is 
invariant under simultanious transformation of the <^'s to ^2 m i (7(i?) mm '^™ where U(R) 
is the unitary representation of SO (3) given by Yi m (Kk) = J2 m , U(R) rnm rYi rn '(k). 

That is, if the BCS equation (1.58) or (1. 62b) of the quadratic mean field model (1.62) 
has a solution such that |A k | 2 is not SO (3) invariant (and, by [FKT2], this is necessarily 
the case for any nonzero solution for d = 3 and I > 2), then the (a ko .a ko -) expectation of the 
quadratic mean field model (1.62) does not coincide with the corresponding expectation of 
the quartic reduced BCS Hamiltonian (1.61). 

In [AB], Anderson and Brinkmann used the quadratic mean field Hamiltonian with 
an £ = 1 interaction to describe the properties of superfluid Helium 3. The basic quantities 
in their analysis are the (a k(J a_ kT ) expectations or the matrix A k(TT which is obtained as a 
solution of the gap equation of the quadratic model. In view of the discussion above, one 
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may regard as the more natural approach to take the quartic BCS Hamiltonian (1.61), to 
add in a symmetry breaking term which breaks the U(l) particle symmetry as well as the 
spatial SO(3) symmetry and then to compute (in the approximation 0™ = 5 q0j o(j) rn ) the 
infinite volume limit followed by the limit symmetry breaking term — > 0. In particular, 
besides the usual U(l) symmetry breaking one may expect SO(3) symmetry breaking in the 
sense that probably also for the (a + a) expectations the above two limits do not commute. 
That is, if B denotes the SO(3) symmetry breaking term, whereas lim^oo limB^o( a k (T Okcr) 
has SO(3) symmetry, hing^o nm L^oo(ak cr «k<r) may have not. However, to what extend 
the quantities of the [AB] paper are related to these expectations is not clear. For example, 
for I > 2 ([AB] has £ = 1, but still uses the quadratic mean field formalism) any nontrivial 
solution of the quadratic model (1.62) is necessarily anisotropic, without any external SO (3) 
symmetry breaking field at all. The quantities liniB^o limL^ 00 (a^ (T ak C r) for the quartic 
model (1.61) or the model discussed in this paper of course would depend on the direction 
of B, if they become anisotropic. A more careful analysis of this question we defer to 
another paper. 
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II. Solution of the Model in the Forward, Exchange 
and BCS Approximation 

Let H = H + H int where 




(1.2) 



and 



l 

L 3d 




(1.3) 



o-,t€{T,I} k,P,q 



Then 




(1.23) 
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where the exponent in the fermionic integral is given by (k = (3L ) 



Vint(V') = ^ Y U ( k -P) A,a^p,a^q-k,rt 

ct,t€{T,I} k,p,q 



q-p,T 



(1.24) 



and k = (k , k) e g(2Z + 1) x M. 

The forward, exchange and BCS approximation is obtained by restricting the above 
sum to the following terms 



forward : Sk, P , exchange : Sk, q - P , BCS : S qj o . 

We consider the approximation Vint ~ Vf orw + V ex + Vbcs = W(V>, ip). That is, 



(1.35) 



W(V>,V0 = i Y U ^ ~ P)[**.P + S k,q-p + Sqfi] ^k,a1pp,a1pq-k,r1pq-p,T (H.l) 

°yr€{t,|} k,p,q 



and 



*7(k-p) = -{ 



1 J 



e^ k e"^ + — if d = 2 



=: -£Wk')^(p') (L37) 
£ ^ A^ m (k')^ m ( P ') if d = 3 z=0 



and we abbreviated 



k = f3L d 



(11.2) 



This model can be solved explicitly. Before we write down the general solution we shortly 
indicate the computation for V nt ~ Vbcs an d U(k — p) = —A. In that case it comes down 
to the usual effective potential computation with 'constant <j>\ That is, using the identity 

((f) = u + iv, 4> = u — iv e C) 



e 2ab = ± J R2 e'Wh-vWdudv 



(11.3) 



one obtains 



Z = J e "(*) H E fe **t*-h x ( a) ii Ep ^_ px 

= / n^detf .; 



feo-ek iX^ (f> 
1 

i\2 (f> —iko — eit 
e -*V (u,v) dudv 



^e-^ 2 dudv 



(II.4) 



R 2 
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In the following theorem the £ variables take care of the exchange contributions, the w 
variable makes the forward contribution quadratic in the fermion fields and the cj> variables, 
as in (11.3,4) above, sum up BCS contributions. 

Theorem II. 1: LetU(ip,^p) be given by (II. 1) and let 

Z(P, L, {s k , a }) = J e^'* + - ^ duett, VO (H.5) 



per I 



K 



pv |s=0 



log Z(/3,L,{s ka }) 



(II.6) 



Let w G R, 
define the fields 



CTT 



^ T + it£ T G C /or < / < J, (ar) G {TT, U, U} and 



j J 

£w = E x hL , s k(JfJ = J2 woo , 

z=o z=o 

1 1 

T kT = t/" 2 w + zS kTT + iS kTT , T ki = t/" 2 w + iH kU + iH kU 



(II.7) 

(II.8) 
(II.9) 



2=0 



z=o 



$w = E toOO - 2/K-k')] , = J2 x i$la MV) - yi(-v)} ■ ( IL1 °) 

1=0 1=0 

where Uq = U(k = 0) = j d d x[/(x). For each k = (ko,k), let S k be the 8 x 8 skew 
symmetric matrix 



Sh = —i 



( 







$kU 





S kU 








c&kTT \ 












$kU 





~ S kU 


$kTT 







-$kU 








~\ A -k[ 





"*k U 


S -kTi 










-$ kU 


\ A -ki 





-$k U 








- S -kU 

















\M 


-^-kn 












5 kU 













-$_ku 


\ 





_$ kTT 







$-kU 










— $kTT 








S -kT4 





$-kU 




/ 


fcer 


= a k - 




= - 


ek - $ka - 


Y ka and let YtS k 


6e t/ie Pfaffian of S 



given by Lemma II. 2 below. Then: 
a) 

Z((l,L,{s k , a }) = ! {(^)VfS fe W(™,£,0) (11.11) 



feo>o 

k<EM 
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where 



dv K (w^,cf>) = {^Ye-^dw n n{(f) 2 e-^^MC^du^dv^da^db 1 ^} 

1=0 (TT K > 

(11.12) 

b) 

d -Of Q 

r si^e-n^,*)^ n <w 

fcW = ? l %,p 7 vi ^77 ( IL13 ) 

where the effective potential V is given by 

= x + EEd^i 2 + i«"i 2 > - ^ E n >n 1^ < IL14 > 

<tt z=o keM 

and nd£# = n n di l aT d$ aT . Here or E {TT, U, U}- 

<7T i = 



Remark: The product n ^ f ^ fc in the effective potential (11.14) where ~PfSk = Pf 'Sk{a,k, a -k) 

k >O a k a -k 

has to be computed according to the rule 

n = lim n wy*"'-*) (n.15, 

fc >0 fc fc e <o fc >0 fc fe 

That there are cases where it is necessary to make the phase factors explicit can be seen 
from the discussion given in the proof of Lemma Al in the appendix. 

Proof: Since we assume U(p — k) = U(\t — p), one has 

$) = ^3 S U ^k,a^p,Ttpk,atpp,T + ^3 ^ ^ ?7(k - p) ^ p ^k,rtpk,atpp,T 

cr,T k,p <j,t k,p 

+ I? S S ^ k ~ P ^ ^P,^-p,Ttpk,atp-k,r 
<j,t k,p 



"^^^(^tV'fcT + ^fciV'fci) (VvtV^t +i>pii>pi / 

k,p 

H ^3- t^(k - p) {^t^PT^r^fcT + ^pl^pl^U^kl + ^p\^pi^ki^k\ } 

h ^ t^(k - p) {$ pT $_ pT V/trV'-fcT + ^pi^-pi^ki^-ki + ^p\^-p\^k\^-k{ } 
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1 v 2 
e 2 X 



We substitute U(k — p) = J2i=o MviO^^Vlip') and use the identities 

-4= Le Xw e-? w2 dw 

-2XY _ 1 r AX^+iYf.- 1 '^' 2 
- 2tt Jt ~ ' 

to obtain 



R 2 



^-2^1 dudv 



ex pfe [< (&) ' 3t £ £ ^(k') ^T^T + < (&) ' ^ T \ £ ^( k ') 
^ i L fc fc 

expfe f z (&) ' *L £ £ woo ^i^i + • (&) ' eh i E ^( k/ ) ^i^i 

^ 2 L fc fc 

expjE [* (£) ' ?u £ E ^( k/ ) + i (£) ' E ^( k/ ) ^i^r 
^ i L /e fc 

ex p{E [* (&) 1 $T i E ^ k ') ^T^-fcT + * (^) ' £ E ^( k ') ^-fcT 
^ i L fc fc 

<*p{£ [< (&) " *U £ E W0O ^-fc| til 71 E ^( k ') M-kl 

^ i L fc fc 

exp{£ f< (£) ' P n £ E ^( k ') + < (£) ' £ E ^( k ') ^rV'-fci 



x 



X 



X 



X 



X 



dv(w, 



where 



= 4=e"^ n n i e-iWUl +1^1 ) du l OT dv l oM ar db l a 



Z (<tt)6 
{TT, U,U} 



By a substitution of variables and collecting terms, one obtains 



-W(V>,V0 



/ ex p{^E ^+<EAf(e} T w(k / )+e{ tW (k / ) 

fc L z 

ex p(^ E k<N + i E A ? (^u ^ k ') + Si ^( k ') 

^ fc L 2 
expj^ z(2A/)^u £ E ^ k ') ^fc| + *(2Aj)'£u J £^(k') ^ fc T 



V'fciV'fcl j x 



x 
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X 



X 



exp ^ i*Mr i E " (k,) T ( " k/) &tVUt + i*M T £ E ^*t^-*t 
i L fe fe 

expj ^ i(2A I )'0' n i^yz(k')^ fcT ^_ fci + z(2A/)V n J J^y/( k ') V'fctV'-fcj. fdzv K (u;,£,0) 

where dv K (w, £, </>) is defined in the statement of the theorem. Using the definition of the 
fields S, T and <E>, the above expression reads 

= J expji Jr kT V'fci^fcT + F k| i^ki^ki + iS kn ^ fcT VfcJ. + ^kTlV'fc^fcT 

We now rewrite the exponent in order to perform the fermionic functional integral. Since, 
if the set of spatial momenta satisfy k e M w = — M^, 

^r fcT ^ fcT ^fcT = E E Ffc T ^tV^t = E E ^T^T^T +r_fc T V'-fcrV'-fct] 

fe fc 6^(2Z+l) k6M„ fc >Ok€M„ 



fc 

fc >o 



one obtains, using the antisymmetry of the &koa, ^kaa 

E{ E Fk<J ^ka^ka + i^U^Pk^-kl + l^U^k^-ki 
k 

+ l^U^k^ki + zS kn V5 fei VfcT + § E^wV'fcaV'-fca + ^aa^ka^-ka) } 
= 5 E 1 E Fk<T i^ka^ka -^kai'ka) + T-ka (^-ka^-ka ~ ^-ka^-ka) 



k 

fc >o 



+ «£ k u (V'fctV'fci - VftiVW) + is_ku (V'-fctV'-fci - V'-fci^-fct) 

+ l E $k<ra(^fc<rV'-fc<r ~ ^-ka^ka) + ^aa(^ka^-ka ~ ^-ka^ka) 



Since, if ajt = iko — 



k,a k,a 
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one obtains 

k,a ak J J k,a 

fe >0 ^ ^ 



dipkidipkidip-k^dip^^ \dv K (w,£, 

7 fe„>0 



ie ; 



where 



and PfSk is the Pfaffian of the 8x8 skew symmetric matrix Sk defined in the statement 
of the theorem. We used that 

II dipkadipka = IT III + dipk^dipk^dip-kidip-kidipkid'ipkid'ip-k'[d'ip-k'[ \ 

k,a fc >0 k ^ J 

Part b) of the theorem follows from 

Lemma II. 2: Let S k be the skew symmetric 8x8 matrix of Theorem II. 1. Then the 
Pfaffian of S k is given by 

+ $kTT^k|4^kn^-kn + ^kn^kU^kU^-kU 
+ ^kTT^kH^kTT^kii 

+ zS kU A_ fcT $ kU $ kU + zS kU A_ fcT $ kU $ kU 

- iEknA-kl^-kU^tt ~ ^ S kT4^-fc4^-kT4^kTT 
+ iH_kt|-Afcj.^ktJ.^kTt + ^--kTi^fci^kTi^kTT 

- iH_kt|-Afct^-kTl^k|| - «S_ kU Afc T $_k U $ kU 

+ SkU S kU^-fcT^-fc.L + S -kT4 S -kT4^fcT^4 
+ SktiS_ kU $ kU $_ kTi + S_ kU S kU $_ kU $kTi 

- S k |jS_ kU $ kTT $ kU - S kU S_ kU $kTT^kH 
+ S kU S kU S_ kU S_ kU 
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Proof: The Pfaffian of Sk is given by the sum of all contractions N-i'P'P) of the fields 



where the value (ipip) is given by the corresponding matrix element. That is, the Pfaffian 
can be evaluated by using Wick's Theorem or integration by parts. Since Sk is an 8 x 8 
matrix, one has Pf[— Sk] = PiS k and 

PfS k =(^ k ^ k] ^- ki 'lp- ki 'lp ki ^ ki 'lp- k ^- k] )s k 

= - A k] {^- kl ^- kl i) kl ^ kl ^- k ^- k} ) Sk - z$kn(V , fcrV , -fciV , fciV'fciV , -feTV'-fcT)s fc 

- ^U^Pk^-kitp-ki^kitp-k^-k]) S k + ^]]{>Pk^-kitp-k{tpki^kitp-k]) S k 

= + A k] - A_ H {ip ki ^p ki ip_ k ^_ k] ) Sk + zA fcT $ kU (ip_ ki ip ki ip_ k ^_ k] ) Sk 
+ iA k] -Z- kn (ip- kl ip kl ip kl ip- k] -}s k 

+ $kU®ku(^kl^kl^-k]^-k]-}s k ~ $kU-ku(^-kl^kl^-k]^-k]-}s k 

- ^ t i^ n (ip- kl ip kl ip k iip-^)s k 

- S kU $ kU (^_ fc ^ A;i ^_ fcT ^_A ;T )5 fe - S kU S kU (^_ fci V-A;iV'-fcT^-A : T)s fe 

- ^ n ^ n -(tp- kl tp- kl tp kl tp- k] -)s k 

+ ^^ n {lp- ki 1p k li) k l1p- k ^) Sk - ^U-^U^-ki^-kilpkilp-k]) S k 

- $ kn $ kn ($_ fci ^_ fci ^ fci ^ fci )s fc 

= + {A k] A- k[ + $ kU $ kU ) (ip ki ^p H ip- k ^- k] ) Sk 
+ (A fcT iS_ kU - ^ n ^ u )(ip_ ki ip ki ip ki tp_ k] ) Sk 
+ (A fcT z$ kU - $ kn Z kn )(ip- kl ip kl ip- k] -i)- k] -}s k 

- Hkti^kTi(^-fci^fciV'-fcT^-fcT)Sfc _ -^\{-^\{ l \i>-k{'P-k{'P-k\ : P-k\)s k 

- ^U^Wi^-klip-kiipkiip-ki) s k + ^kn^kU^-fc^fc^fciV' — k]'/S k 

- i kTT H kTi (^_ fci V-fciV'fciV'-fcT)Sfc ~ ^kTT $ kTT(^-H' i / , -fciV'fc^fc4)5 fe 
= + {A k] A- k[ + $ kU $ kU ) {A ki A_ k] + $_ kU $_ kU ) 

+ (A fcT iS_ kU - <5 kn $ kTT ) (-$ kU $_ kU - iE_ kn A ki ) 

+ (A fcT i$ kU - $ kU S kU ) (-z$ kU A_ fcT - E_ fcn $_ kn J 



- S k |j$ k |j (-z<& k .u^4-fcT - H_ktj,^-ktl) 

- S kn H kn (-A_ kl A_ k] - - S_ kU S_ kU ) 

- Skti^ktT (^-fej.*^-kn + $kU S -fcu) 
+ ^kTT^kTl ($k±j$-kU + lE-ku^fcJ.) 
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- 5>kTT $ kTT ( - A- H A H - $ kU $ kU ) 
By multiplying out the brackets one obtains the stated formula • 

Before we specialize Theorem II. 1 in section III where the effective potential and 
the two point functions are computed more explicitly, in the following theorem we write 
down the integral representation for the generating functional of the connected amputated 
Greens functions. 

Theorem II. 3: Let U '(ip , V>) be given by (II. 1) and let V be the effective potential (II. 14). 
Let 

G(rj) =log J e- w( ^ + ^ + ^ c (V>,V0 (11.16) 

be the generating functional for the connected amputated Greens functions. Then, if C k 
denotes the eight component vector 

C k = (akVkh -a-kVkha-kV-ki, -a—kV-ki, -a—kV-kh a -kV-kf) 

and Sk = Sk(w,£,4>) is the 8x8 matrix of Theorem III, one has the following integral 
representation 

G(q) - G O = -\ > a kVka r] ka + log ^ — — ^ 11.17 



Proof: By a substitution of Grassmann variables, 

e -u(i,+ v ^+fj) d/Jc = 



e ~ « S fcCT a kVk<rVk<r / e -W(V>,"0) pj _k_ e T, J2k* «fc[V'fc<T'7fc CT +'7fc CT ^fc<T-V'fc CT ^fe<T] jj dtfjkadlpka 
J ka ah ka 

As in the proof of Theorem II. 1, one has 
such that 

-U(i>,i>) _k_ e \ J2 k<7 a k[4>kaVk«+Vk«4>k«-4>k«->Pk<?} dlfjkadlpka 
ka ak ka 
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where 




Yl ( Pfgfc ^ e \ Z]fe CT a kVkkaVka+Vka^k^] ^ 

fc >0,k V a fc a -fc/ 



//,> n o, k (^;) eiEw< "' £ * >x 

fc >0,k ri,5fc fc >0,k 

y fc >o,k \ k ~k / 



By definition of F 

>0.k V a fe a -fc / 



fc >0,k 

which proves the theorem • 



III. Solution for Pure BCS 

We consider the model 

Z(P, L, {s k , a }) = J e - U( ^ )+1 « Z>,. s ^^ d ^ $ (IH.1) 

with e k = k 2 /2m — fi, C(k) = l/(iko — e k ) and 

U{#,$) = ^ Y. U ^-^)^,^-k^ P ,^- P ,r (HI.2) 

The electron-electron interaction U is given by (1.37). 

To write down the effective potential and the two point functions in this case, one 
first has to compute the Pfaffian of the matrix S k of Theorem II. 1. Since we consider only 
a BCS interaction, the S and V fields are zero. With Lemma II. 2 one obtains 

PfS fc ={a k] a- ki + $ kU $ kU ) (a- k] a kl + $_ kU $_ kn ) 
+ a fcT a_ fcT $ kU $ kU + a fei a_ fci $ kn <i kTT 

+ $ kT T$kU$kU^-kU + ^kn^kU$kU$-kU + $kn$kU$kn$kU ( IIL3 ) 
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where a ka = &k ~ s ka- In particular, for s ka = 

PfS k = (a k a- k + n+)(a k a- k + Q") (III.4) 
where are the solutions of the quadratic equation 

O 2 - ($ kn $ kn + $_ kn $_ kU + $ kTT $ kTT + $ kU $ kU )o + $ kTi $ kTi $_ kU $_ kTi 
+ $ kn $kU$kU$-kU + $kTT^k|4$kT4*-kT4 + $kTT$kU$kn$kU = (III.5) 



The effective potential becomes 

TV 

fc o >0 



"(0 = E E i^i 2 - ^ E "g.n } 

<"-€{Tt,U,U} 1=0 k6M 



We start in section III. 1 with the case of a delta function interaction, that is, J = in 
(1.37). This reproduces the usual mean field results. In section III. 2, we consider a more 
general electron-electron interaction of the form (1.37) with arbitrary J. One finds that 
the standard approach based on approximating the Hamiltonian by quadratic terms and 
imposing a self consistency equation may be misleading. 

III.l BCS with Delta Function Interaction 

Corollary III.l: Let s k ^, Sk,i, r k , r k be some real or complex numbers, let 

Z(p, L, {r k }) = / exp<^ £ ip k ^_ ki ip p ^- pi + ± H[s k ^ k ^ k ^ 

J I k,p k 

+ Ski^ki^Pki + r k ^k\^-k{ - fk4>k^-k i] ^d/j,c(ip, (HI.7) 

and let 

i(Vw/vW = ^logZ(P,L,{ 8k },{rk})\ ah=0frk=r (IIL8) 
^ P ^-Pi)p,L,r = ^logZ(/3,L,K},{r fc })| Sfc=0rfe=r (III.9) 

Define the effective potential 

Vfc r («, v) = + (v + M)2 -f M *± Wp M (\(u* + v*)) (111.10) 

where 

W^(y) = ±\og 



B n 2 

cosh (? y e k+^) 



cosh~|~ek 



(in.ii) 



25 



a) Let a p = ipo — e p . There are the two dimensional integral representations 

Z(f3, L, r) = * J e~ KV ^^ u ^dudv (111.12) 



_/■ 






v ">dudv 


a p a- p -\-\(u 2 -\-v 2 ) 






v ">dudv 




_/ 


-iV\e' ia (u-iv) 


— K,Vf3 t r(u 


' v ^dudv 


a v a- v -\-\(u 2 -\-v' 2 ) 






' V ^dudv 




_/ 


-iV\ e ia (u+iv) 




' v ^dudv 


a p a- p +\(u 2 +v 2 ) 






• V ^dudv 





WM^, - (ni.i3) 



Uj> v t^_ p I )B L r = a * a -*+^+ v ;> f (IH.14) 

k \^p,T^ P,i/^,r f e - KV 0,r(^)dudv y 1 

b) Let V£(p) := p 2 -f -(0dWp M (\p 2 ) and |A| 2 = \p 2 be given by Vp(p Q ) = min p > V f3 {p) 
Then, if A = A (A, (3) = \ A| e~ ia , one has 

lim lim ^logZ(p,L,r) = -V P (\A\) (111.15) 

iJ^oi^o ih^v^P,*)^ = ap ~_ a p +\ A \ 2 ( IIL16 ) 

iJj^oiSo JU^P^-P^P^r = 0pQ _ p A + | A p ( IIL17 ) 

l^oi™, ^$P^-P^P> L ' r = a p a~t\A\ 2 ( IIL18 ) 

lim lim J(V'p,tV'-pa)^,r = r lim }\ m n i(Tp P ,^-P,l)p,L,r = (111.19) 



Proof: a) For a 5 function interaction one has J = in Theorem II. 1 and the fields (11.10) 
are zero, $ kfJ(J = 0. The fields (II.9) become $ kn = (2\) 2 (p° tl = (2A)*0. Since we added 
the and fiptfj terms to the exponent in (III. 7) which were not present in Theorem II. 1, 
these fields have to be substituted by (g = V2X) 

g<f>-> gcj>- ir k , g(p^ gcp + if k 

The Pfaffian (III.3) becomes 

Pf S k = ( a fcT a_ fci + \g<j> + if k ] {g<j> - ir k ] ) ( a_ fcT a fci + \g<j> + if _ fc ] \g<f> - zr_ fc ] ) (111.20) 
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such that 



r?— PfS. 



V 



PfS j, | Sfe=0irfc=r a p a_ p + [g<f>-ir][gct>+ir] 

and 

~ Pf,g P _ -i[g4>+if] 



Pf5 P | Sfc =0irfc =r a p a„ p + [g4>-ir][g<f)+ir] 

The effective potential (III. 6) is given by (the \(p\ 2 term is not shifted by r) 

I aka-k + [g<t>-^r\ [g<p+ir\ { " 

fe o >0 



(111.21) 
(111.22) 



V(<t>) = \4>\ 2 - ^ £ lQ S n { ^±^glfeM } 2 (111.23) 



The product over /c € ^(2N + 1) in (III. 23) can be computed explicitly using the formula 
[Hn] 

oo / , \ cosh f f v/oHT) 

n o (i + (4^) = c l o 2 sh fa j (ni.24) 



One obtains 

fc n {^?} = fco6f g N+1) ( 1 + 4^) = { COS c h osh(tSf } } ^ IIL25 ) 

which gives, if one approximates the Riemannian sum X^keM^ by an integral 

= \<P\ 2 - [ ^W p , k ([gcf>-ir}[g4> + if}) (111.26) 

Thus one arrives at the integral representations 

Z(0, L, r) = * J e- KV ^dudv (111.27) 

^p^p,ah,L,r = fe-^Wdudv (IIL28) 

f -i[gt+ir] e ~ KV ^dudv 

1 /„/, „/, \ J apa- P + {94>-ir]{g<t>+tr] /ttt 00 n 

« ^- P ,l)?,L,r = f e -«V Wdudv ( IIL29 ) 

where V is given by (III. 26). Part (a) then follows from the substitution of variables 

J F(g(f)-ir,g4> + if) e~ K ^ dudv = J F (e lQ ^, e~ %a g$) e ~ K ^ HvJr ^ f ) dudv (111.30) 

which holds for both signs of A. 

To obtain part (b), one has to compute the limit of 

e —KVp, r (u,v) 

^,rM= Je _ KVprMdudv (IIL31) 
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where 



Vfc P («, = u " + (v + ^) 2 - W P {X(u 2 + v 2 )) 



(111.32) 



if we define 



M 



M 



cosh(-^e k ) 



(111.33) 



In particular, 



signW / 3(A('U 2 + n 2 )) = signA 



For positive A and nonzero r, Vp tr (u, v) is real and has a unique global minimum determined 
by 

2u-2uXW' p (X(u 2 + v 2 )) = 
2(v + \r\) - 2vXW p (X(u 2 + v 2 )) = 

Since XW' p = 1 does not solve the second equation, the only solution of the first equation 
is u = and one is left with 

v [XWp(Xv 2 ) - 1] = \r\ 

which has a solution v = 0(\r\) which is a local maximum and two nontrivial solutions 
Af 2 = |A| 2 + (9(|r|) where the positive one is only a local minimum and the negative one, 
vq, is the global minimum. Therefore 

e —K[V0,r(u,v) — Vp, r (O,V O )] 

lim 5 K r (u, v) = lim -z tt-. — — r— r . — v-- — ^ = S(u) S(v — vo) 

k^oo ' v ' ' K ^oo J e - K l v pA u > v )- v i3A > v o)]dudv 



and 



lim lim S K r (u, v) = 5{u) 5 { v + — ) 

|r|— >0«— >oo ' \ 9 / 

This proves the formulae under (b) for attractive A. Since 

lim V rj p(u, v) = u 2 + v 2 - Wp (X(u 2 + v 2 )) 



(111.34) 



M— 



is an even function in u and v, lim r ^ (' i / ; p,T' i / ; -p,J.)/3,L,r=li m r^o(' i / ; p,T' i / ; -p,4)/?,-L,r = 0. The 
limit of the logarithm of the partition function becomes 

i logZ = i log * I e- KV e^ u ' v) dudv 

= £log* J e- K[v ^ r{u ' v) - v ^ Ao ' Vo)] dudv + ±\oge- KV e- r(0 > Vo) 
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d 2 V 

The first term on the right hand side may be approximated by (V uv = du ^ v r (0, vq) = 0) 
ilog^ [ e-^ v ™ u2+v ^ v - v ^dudv=±\og± [ e-^ v - u2+v ^dudv K ^° 

re & 7T / re & tt / 

which results in 

lim llogZ = -^, r (0,t;o) 

Now let A be negative. In that case the effective potential (III. 32) is complex: 

V Ptr (u, v)=u 2 + v 2 -^- 2iv-$- - Wp {\(u 2 + v 2 )) 
= p 2 + \W,(\p 2 )\-^-2ipcos^ 

Since the real part Up jr = ReVg jr has a global minimum at u = v = one has, since 
£/"=^(0,0) = ^(0,0)>0, 



-re [l/ /3ir (« I t;)-2it;-Ji.] -re (uW)-2™^Ll 



Je ^ } dudv Je 12 ^ } dudv 



U" 2 -Ju^ v 2_ 2iv JfL] 



e -re^-u e 2 L ^T|C7»J 



which results in 



T777 F - ► O(tt) f — Z H= 



lim lim 5 K r (u, v) = S(u) S(v) (111.35) 

| r U0 k^oo ' 



and 

|i}™ J™o i($P^-P,l)(3,L,r = J apa _~+x&+v*) S(u)6(v) dudv = ^ = 
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Using (11.17) and (111.34,35) the infinite volume limit for the generating functional for 
the connected Greens functions can be computed in a similar way. One finds 



Corollary III. 2: Let r = \r\e ia , 7 = ge ta (u + iv ), 7 = ge ia (u — iv) and let 

G(n) = G r (P,L, v ) = log J e w ^+^+^c(V^) (HI.36) 
be the generating functional for the connected amputated Greens functions where 

U r (ip, $) = -^s^2 ^k^-ki^p^-pi + £ S ^[rip k] ip- ki - f^fc T ^_fcJ (111.37) 

k,p k 

Let Vp !r be the effective potential (III. 10) and let A be given as in Corollary III.l. Then 
there is the integral representation 

G{rj) - G(0) = - J [ a fc%T%T + a -kV-kiV-ki\ + (111.38) 
k 

+ loe — 

J dudv e - KV eAu,v) 

For attractive A > the infinite volume limit of the generating functional is given by 

lim lim {G(P,L, r, rj) - G(fi, L, r, 0)} = (111.39) 

\r\^0 L^oo 

1 T |A| 2 - |A| 2 

- k 2s [°* a k a J+\A\- + O-fc afca l fc V|A|^ V-HV-kl 

k 

+ A a fc a_ fc + |A|* V-HVkt + A afca _ fc + | A |. V^V-kl 



For repulsive A < one obtains 



lim lim {G(0, L, r, 77) - G(f3, L, r, 0)} = . (111.40) 

H— >0 L^oo 



III. 2 BCS with Higher Angular Momentum Terms 

We now consider the case where the electron-electron interaction contains higher an- 
gular momentum terms. In this case one finds that the usual mean field approach, based 
on approximating the Hamiltonian by quadratic terms, may be misleading. To simplify 
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the algebra we assume that only even £ terms contribute in (1.37). In that case the fields 
^kcro- with equal spin (11.10) are still zero. So let 



U(k'-p') = { 



£ =■■ E *viWvi(p') ( IIL41 ) 

E E XeYim (k') Ytm (p') if d = 3 i™„ 



fco m =- 

V. £ even 



where y^( — k) = (— l)'j//(k). We only consider the expectations, so we let r = 0. The 
Pfaffian (III. 3) becomes 

PfS fc = (a fcT a_ fci + $ kn $ kn ) (a_ fcT a fci + $_ kU $_ kn ) (111.42) 

where 

J J 
d> kU = £ (2A,)»0' n wCk 7 ) , *ku = E (2A ; )^ Z Ti ^(k') (IH.43) 

1=0 1=0 
I even / even 

and the effective potential is given by, using (III. 25) again to compute the ko product 

J . 

V(<hl) = E I^Tll 2 " / (Cf W ?A^U^U) (IH.44) 
i=o Jm - 

where Wp^ is given by (III. 11). By Theorem II. 1, the two point function is given by 

±(ipp*ippv) = j (111.45) 

/ e- Kl/ (^Tx) n du l n dv l n 

To compute the infinite volume limit, one has to find the global minimum of the real 
part of V(4>ii). This is easier for d = 2. However, using symmetry arguments, it is possible 
to give a rather explicit expression also for d = 3 without knowing the exact location of 
the global minimum. 

Consider first the two dimensional case. An analysis done by Albrecht Schuette in his 
Diploma thesis [Sch] shows the following result: 

Suppose that A m > is attractive and A m > for all £ ^ m. Then 

lim ^t/wVwW = ctxI 2 2 (111.46) 
where p m is determined by the BCS equation 

rf|k||k|^^^ = (111.47) 
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The form of the two point function (III. 46) can still be obtained by applying the stan- 
dard mean field formalism [AB,BW]. However, the situation is different in 3 dimensions. 
Before we state the corresponding theorem, we shortly recall the mean field equations 
(1.57,58) 

The (a + a) expectations are given by 



„ + n \ _ 1 / 1 p tanhd^+A^Ak) ' \ , , 

a ka°ka) - 2 I 1 - e k ,„ 2 , A . A . ) U-5') 



where the 2x2 matrix A k , A k = — A_ k , is a solution of the gap equation 



f dgk U{ , _ k /) A tanh(| ye k+ A- AQ 

JM 1 ^ 2v /e k +A*A k 1 > 



In 2 dimensions, if one substitutes U(p — k) by a single attractive term A^e^^ p 
then (1.58) has the unitary isotropic solution A k = A ( _ e « Vk J for even t and A k = 

A ( cos ^ k sin ^ k ) if £ is odd. In that case A* A k = \A\ 2 Id and (1.57) coincides with 

V smtipk — cos£<^ k J k K l l V / 

(III. 46) (after integration of the latter over po)- 

In 3 dimensions, it is proven [FKT2] that for all £ > 2 (1.58) does not have unitary 
isotropic (A k A k = const Id) solutions. In view of that result, the symmetry considerations 
below indicate that in 3 dimensions for £ > 2 the standard mean field approach is misleading 
since one would no longer expect 50(3) invariance for the (a k(J ak(j) expectations according 
to (1.57). But this is indeed the case if there is no external SO (3) symmetry breaking term 
(which is also not present in the quadratic mean field model, see the discussion following 
(1.61) at the end of the introduction). 



Theorem III. 3: Let £ be even, Xg > be attractive and let 

i 

U{$^) = j±l_Y. E *Wk')*Wp') M-kiM- P i (HI.48) 

k,p m= — i 

Then, if e# k = e k for all R G 50(3), one has 

lim I(VS p ^ pCT ) /3jL = lim / \^ pu e- u ^^dii c {^^) 

= I iPQ+ej, df2(x) (m 4q x 

J s2 P g+e|+A^jj|S m aO i y /m (x)P 4tt 

where po > and a E C 2i+1 , £ m |a^J 2 = 1, are values at the global minimum (which is 
degenerated) of 

V(p,a)=p 2 - [ (0ilog 

JM 

32 



cosh(JVi|+A^gj S m Q m y £m (k')P) 
cosh 



In particular, the momentum distribution n p is given by 

n » = jfe ^W^/^ = J s2 2 " e P T^aCT J T (IH.50) 
and /ios 50(3) symmetry. Here A(x) = A|po S m o^^mC*)- 

Proof: Substituting tt^ by u m and by — v m in (III. 45), one has to compute the infinite 
volume limit of 



p -KV(</>) yi£ J ? . J., 

c ll rrl = _/) uu> m uu r 



/ R «+a e n^ = _^ du m dv m 



(111.51) 



where 



and 



n*)= E i<u 2 -/ (0f log 



cosh 4 e k 



$k = A| ]T m ^ m (k') 
Let C/(i?) be the unitary representation of 50(3) given by 

mm' m' \^ ) 

m' 

and let E m (W)m*W k ') =: (L/$) k . Then for all i? G 50(3) one has (C/(R)$) k = $ fl -i k 
and 

cosh ^e k 



K(£W) = £ |[l/(fl)# m | 2 - y jf ^ § log 



= X>-i 2 -/ M (fip? lo « 



cosh 



cos 



cosh~^~e k 



= 1/(0) (111.52) 



Let 5 4 ^ +1 = {0 G C 2e+1 | Em \0rnl 2 = 1} ■ Since C7(R) leaves S 4e+1 invariant, S 4i+1 can 
be written as the union of disjoint orbits, 

where [a] = {U(R)a | R G 50(3)} is the orbit of a G S 4£+1 under the action of U(R) and 
O is the set of all orbits. If one chooses a fixed representant a in each orbit [a], that is, if 
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one chooses a fixed section a : O — > 5' 4£+1 , [a] — > <7[ a ] with [<7[ a ]] = [a], every e C 2i+1 
can be uniquely written as 

cj> = pU(R)a [a] , p=\\(f>\\>0, <x=^, [a] eO and ReSO(3)/I [ce] 

where 7[ Q ] = = {S 1 G 50(3) | L r (5')(j[ Q ] = cr[ a ]} is the isotropy subgroup of <7[ a j. Let 

/ R «+a Tim dumdvm f((j>) = J R+ Dp J Q D[a] J [a] DR f (p U(R)a [a] ) 

be the integral in (111.51) over R 4 ^" 1 " 2 in the new coordinates. That is, for example, Dp = 
p A( - +1 dp. In the new coordinates 

2 I 2 

\® P \ 2 = \ep 2 \Ys( u ( R )°n) *Wp') = a, p 2 ^ [a] , m r, m (irV) 

I — \ / m I — 

m m 

such that 

a p a_ p + |$ p | 2 = a p a_ p +A^ 2 |E m <T [ J, m y £m (i?- 1 p')| 2 = / (z 9 ' W ' ^ *P) 

Since V(0) = V^(p, [a]) is independent of i?, one obtains 

f ~ a ZL i2 e-" v W n du m dv m 

-(VvVV/ = ^ 

/ e-" v M n du m dv m 

m= — £ 

J R+ Up / Q D[a] J [a] DR f(p, [a],R-ip) e~^M) 

/R + ^/o^[«]/[a]^^^' [a]) 
/ R+ ^J D[a] vol([a]) J -& j . ^ e -*v(p,[«]) 

= / R+ Dp/ £>[a]vol([ai)e-^(p.W) (IIL53) 

It is plausible to assume that at the global minimum of V(p, [a]) p is uniquely determined, 
say po- Let O m i n C be the set of all orbits at which V(po, [a]) takes its global minimum. 
Then in the infinite volume limit (III. 53) becomes 

f DR f(p,[a],R- 1 p) 

L D [a] vol( [a] ) ^ 

hm ±<VvVv> = f nr i m i"' (HI.54) 

k^oo J G . -D[aJ vol([aJ) 

Consider the quotient of integrals in the numerator of (III. 54). Since 
f(p, [alR-'p) = f(p 2 \Xv [a] , m Ytm (^"'P')! 2 ) 

\ m / 



= f[p 2 \T,{U(S)a [a] ) m Y em {R-'p^l 
= f(p 2 \Za [a] , m Y eni ((RS)- 1 ^ 2 ) = f(p, [a], (^)" 1 p) 

34 



for all S G J[ a ], one has, since [a] ~ 50 (3) //[ Q ] 

/ I2R /(p, [a], R-'p) Iso(3)/i [a] DR /(P. H ^P) /j [a] ^ 



/fal ^ 



/sO(3)/J [a] ^//[a] ^ 
/SQ(3)/J M ffl J/ M 13 g f(P> H ( i? ' g )" 1 P) 
IsO(3)/I [oi] Il la] DRDS 

Iso(s) DR f^ t«]^ _1 P) 



/SO(3) x _ p 



/SO(3) 

J 52 dfl(x) | sof3 ,_„ £W f( P , [a], iT x p) 



DR 



J s2 rfO(x) /(/?, [g],x) J go(3) _ p 
J S 2 dfi(x) / so(3) _ p £W 



(111.55) 



where 5 , 0(3) x ^ p = {.R G 50(3) | i?x = p}. If one assumes that DR has the usual 
invariance properties of the Haar measure, then f so ^ DR does not depend on x such 
that it cancells out in (III. 55). Then (III. 54) gives 



L. D a vol( a ) ^ — r — 

J o min l J U \) f dQ(x) 
lim -(ippalppa) = f p-j^ r-fy 

k-oo« J 0min D[a] vol([a]) 

Now, since the effective potential, which is constant on O m ; n , may be written as 



(111.56) 



V(p, [<*])=[ ^0(p 2 |E m a [Q])m y £m (x) 
Js 2 v 



with G(X) = p 2 - J f£ log 



cosh ^e fe 



, it is plausible to assume that also 



J s2 rfO(x) /(p, [q],x) 
J S 2 dft(x) 



_ /" dn(x) 



jgO+gp 



e p +A<.p 2 | ^ m (x)| : 



is constant on O min . In that case also the integrals over O m i n in (III. 56) cancel out and 
the theorem is proven • 
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IV. Solution with Delta Function Interaction and 
Its Perturbation Theory 



IV. 1 Solution with Delta Function 

We consider the model 



Z(P, L, {s k , a }) = J e -"<**> + * °^^ d ^ $ (iv.l) 



where 



U(lp,1p) = ^ + *fc.9-P + 5, ?,o] ^ka^q-kr^pa^q-pr 

= A + + S *>°1 ^k^q-kilpp^q-pd (IV.2) 

It is useful to explicitly cancel the Tpfipfipfip^ and V'j.V'j.V'J.V'J. terms before directly applying 
Theorem II. 1 because then the integral representation becomes 4 dimensional instead of 5 
dimensional. Both representations are of course equivalent. One obtains 

9 pfff 

/ "pis " e-" v ^''"^dv dw xdx f>df> 

(tMw = /Ji ^, — J ;-.v,(^ J . )(fa , (to;nfaM , — < IV3 > 

where z = v + iw, 

PfS k = (A^A- kl + \p 2 )(A- k] A ki + Xp 2 ) + (A k ^A H + A- k ^A- ki )\x 2 

+ 2X 2 x 2 p 2 + X 2 x 4 (IV.4) 

A k ^ = a k - s fcT - y/Xz, A kl = a k - s kl - \f\~z 
as before a k = iko — and the effective potential is given by 

v (3 (z,x,p) = \z\ 2 + x 2 + p 2 - 1 ± s J2 lo s n ffi- ( IV - 5 ) 

The product over ko is computed in Lemma Al in the appendix. One finds 

V (3 (v,w,x,p) = (IV.6) 

" cosh 2 (f V(e k +VA,)HA P 2 ) - sin 2 (| v /A(w 2 +x 2 )) _ /? ^ 



f 2 + «; 2 + z 2 + p 2 - ^ ^ log 

k6M„ 



cosh 2 -2e k 



The 2 variable sums up forward contributions, a; = |£| sums up exchange contributions 
and p = \(p\ collects the BCS contributions. Pure BCS is given by z = x = and (IV.6) 
coincides with (1.47) or (III. 10) (for r = 0). 
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To compute the infinite volume limit of (IV. 3) one has to find the global minimum of 
the effective potential (IV. 6). 

Consider first the case of an attractive coupling A = g 2 > 0. To make (IV. 6) small, 
the numerator in the logarithm has to be large. Hence for A positive w 2 + x 2 has to be 
zero. This gives 

Pf.<7_ 

a-p—gv 



b p1 



PfS p {a p -gv][a- p -gv]+g 2 p 2 (IV.7) 

As in section III.l p is positive for sufficiently small T = 4 and v may or may not be 
nonzero to renormalize e p . Thus for attractive coupling the two point function becomes 



lim ^(V'ptV'pt) 



ip +e pg 



(IV.8) 



if e pg = e p + gvo, \A\ 2 = Ap 2 , and vo, po are the (A dependent) values where V takes its 
global minimum. 

Now consider a repulsive coupling A = — g 2 < 0. In that case the numerator in the 
logarithm in (IV. 6) reads cosh 2 (|^/(e k + igv) 2 — g 2 p 2 ) + sinh 2 (^g^/w 2 + x 2 ) and the 
value p = is favourable. Furthermore v = seems favourable. In that case 



-PfS„ 



a p — gw 



a p— g 2 (w 2 +x 2 ) 



(IV.9) 



and the two point function becomes, if 5 2 = g 2 (w 2 + x 2 ) > 



lim ±(VvtVm) 



a%-6* ~ 2 



+ 



ipo—e p —Sg ' ip -e p +5 



(IV.10) 



since the effective potential is even in w and therefore the w term in (IV.9) cancels. At 
zero temperature, this results in a momentum distribution 



Jim ^(a+a k )f3 iL = lim / f 



dkp 1 
ir 2 



+ 



e ^ t(ik -e k )-8 g 1 e lfc 0^(ifc -e k )+5 g _ 



X(e k - 5 g < 0) + x(e k + 5 g < 0)] 



1 if e k < —S g 
\ if -S g < e k < 5g 
if (5 9 < e k 



(IV.ll) 



Thus, if 5 g is nonzero, the exchange and forward contributions lead to a splitting of the 
Fermi surface. A nonzero 5 g can be achieved by making the coupling g 2 sufficiently big: 



Let K = l 4^J-mk^~ 2 ~ k d F /p and let 



(27r)° 



Vp(y) = y 2 



/ {2^Y S3 



log 



cosh 2 ( | e k ) + sinh 2 ( § gjy) 
cosh 2 4 e k 



(IV. 12) 
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Then 



lim V p (y) = (1 - Kg 2 )y 2 + K[gy - u]\ (IV.13) 

p^oo 



where [x] + = x for positive x and otherwise. The global minimum of is at y m ; n where 

v ■ -1° if 2 2 <V^ (IV 14) 

ymm "\Kc^7 if^7 2 >l/K 1 J 



IV. 2 Perturbation Theory 

In this section we consider the perturbation theory of the model 

Z ex/BCS (P,L) = J e-^/BcsW^)^^) (IV .15) 

For BCS, one obtains a power series in C(k)C(—k) and for V ex one gets a power series in 

C{k)\ 

We will find that the linked cluster theorem, log Z is given by the sum of all connected 
diagrams, cannot be applied to (IV.15). More specifically, whereas the series for Z(X), the 
sum of all diagrams, converges for sufficiently small A (for finite (3 and L) with Z(0) = 1 
which implies that also log Z is analytic for A sufficiently small, the sum of all connected 
diagrams has radius of convergence zero (for finite f3 and L). 

We start with BCS. First we show how the integral representation (III. 12) (for r = 0) 



r - K (/-^E k lo s 
Zbcs = 2k e 



TP 

cosh ej^ 



pdp (IV.16) 



is obtained by direct summation of the diagrams without making a Hubbard Stratonovich 
transformation. That is, without using the identity 

•\ \- r . . . . In- . ,2 



One has 

Zbcs 



= J e-^^ M - klM - pl dtic(M) (IV.17) 

oo / A \ n 

= E^r E ^[^*i^)]i<ij<„<^[^*i^(-**)]i<*j<„ 

n=0 *l,"',*n 
PI,--- ,Pn 

OO 

= E(i) n E e - E C{k 1 )C{-k 1 )---C{k n )C{-k n )5 kuKl .--5 kn , Kn 

n=Q ir€S n ki,---,k n 
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This is the expansion into Feynman diagrams. It can be summed up by collecting the 
fermion loops: 

Say that the permutation n is of type t(n) = l bl ■■■n bn if the decomposition into 
disjoint cycles contains b r r-cycles for 1 < r < n. Necessarily one has lb± + ■ — h nb n = n. 
The number of permutations which have b r r-cycles for 1 < r < n is 



61! • • -bj l bl • • - n 6 ™ 
The sign of such a permutation is given by 

t ^ _ (_ 1 ^(l-l)6 1 + (2-l)6 2 + --- + (n-l)6„ _ {_Xf~Y2=x K 

Therefore one obtains 

00 n . 

nl 



z = E(D" E 



n=0 



61! ■■■b n \ l 6 i • • -n b " 



(_i)"-Er=i^ x 



IblH Ynb n =n 



= E E "'nnHER^^r 



n=0 bi,"-,t„=0 r=l 

lbl H Ynb n =n 



(IV.18) 



The only factor which prevents us from an explicit summation of the above series is the 
n\. Therefore we substitute 



n\= I 
Jo 



e x x n dx 



and obtain 
Z 



^E E nnBEI-^cwct-*)] 



n=0 bx,--,b n =0 r=l 

lblH Vnb n = n 



dx 



-k) 



dx 



poo 00 op 1 C 

" / ° r=l b r =0 ^ 1 fc 

= r^-E.EZi^-^^w 



V log 



l+^C(k)C(-k) 



2k 



2k 



- dx 

) 



o +e kJ/ p<ip 

cosh(^^/ e 2 +Ap2) 
cosh ^ e k 



(IV.19) 
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which coincides with (IV. 16). In the last line we used (III. 25) again. The case of an 
exchange interaction is treated in the same way: 



E C{k 1 ?-..C{k n ?8 kl , Kl ...8 kn , Kr 

TtZiSn ki ,k n 

b r 



n=0 

oo 



= E E »>II^B£Kcm 2 ] 



n=0 b 1: --,b n =o r—1 

lbi H \-nb n =n 



poo 

= e 
Jo 

r 00 -k 
= 2k I e 
Jo 



k Z_-/r=l r 



-^X C(k) 2 



y -jd 



cosh 2 ( ^ e k ) — sin 2 ( ^ \/A y) 
cosh 2 lr e k 



(IV.20) 



(IV.21) 



where in the last line Lemma Al in the appendix has been used to compute the ko product. 
Observe that lim^^ ± £ fc C(k) 2r = for all r > 1 since / ffi (ifco l_ e) j = for all j > 2, 
but for large repulsive coupling linig^oo logZ ex = -V^l/min) > by (IV.13,14). 

We now consider the linked cluster theorem. It states that, if the partition function 



z w = E An E val ( G ) 

n=o Ger n 



(IV.22) 



is given by a sum of diagrams, T n being the set of all n'th order diagrams, then the 
logarithm 



logZ(A) = ^A" val ( G ) 

n =o Gere, 



(IV.23) 



is given by the sum of all connected diagrams. This theorem is easily illustrated for a 
quadratic perturbation. Namely, if 



EH^ E det^S^Cih)] 

l<i,j<r 

it' 



fel • ■ • fen 



= E E nnl-S»i 



n=0 b 1 ,--,b n =0 r —l 
lb 1 H hrtb„ =" 

e E fe log[l+AC(fc)] 



(IV.24) 
(IV.25) 
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then the sum of all connected diagrams is obtained by summing all the terms with b\ = 
■ ■ ■ = b n -i = and b n = 1 in (IV. 24). That is, one gets 
00 i r i 1 

£ 1! H £ l- XC ( k T = £ ^g [1 + XC(k)\ (IV.26) 

n=l ^ k k 

which coincides with logZ. Now consider Z-qcs- The sum of all connected diagrams is 
given by the sum of all the terms with b\ = ■ • • = b n -i = and b n = 1 in (IV. 18). That 
is, one obtains 

oo 

n=0 k 

which has radius of convergence zero, at finite temperature and finite volume. How- 
ever, Zbcs(A) has positive radius of convergence and Zbcs(O) = 1 which means that also 
logics (A) is analytic for sufficiently small (volume and temperature dependent) A. That 
is, in this case the linked cluster theorem does not apply since the right hand side of (IV.23) 
is infinite. 

We remark that we think that this is an artefact of the specific model at hand for the 
following reason. Suppose for the moment that the k sums in (IV. 17, 18) are finite with 
iV different values of k. Then Zbcs(A) is a polynomial in A of degree (at most) N. In 
particular, the coefficients 

£ Iln H£ [c(A0C,( ~* )3 r =0 ifn>iV - ( IV - 28 ) 

!>!,■ ■-,b n =0 r=l r ' ^ k ' 

lf>lH ^nb n =n 

That is, there are strong cancellations between fermion loops of different orders. However, 
for Vi n t ~ Vbcs (or forward or exchange), an n'th order connected diagram is necessarily a 
single n'th order fermion loop such that there are no cancellations at all. In a more realistic 
model a connected diagram contains fermion loops of different orders and cancellations 
are present. In fact, a careful diagrammatic analysis [FT,FKLT,FST] shows that the 
renormalized sum of all connected diagrams of the two dimensional electron system with 
anisotropic dispersion relation (such that Cooper pairs are suppressed) has positive 
radius of convergence (which is, in particular, independent of volume and temperature). 



Appendix 

Lemma A.l: a) Let = ik — e^, e_k = and let b,c,d be some complex numbers. 
Then 

n a k a- k +ba k +ca_ k +d ._ j. TT a k a_ k +be lek o a k +ce~ lek ° a_ k +d 

«fc«-fc ' e^O J-J- «fe«-fe 

fcoef (2Z+1) *<o fc 6f (2Z+1) 
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_ cosh 1^+0;+ coshf v /e2+^_ ^ {b+c) ^ ^ 

cosh ^e k cosh |Je k 

u>/iere u± = p± Vp ~ 4 ? w ahp = (6+c-e k ) 2 -e^-46c+2<i and g = 46ce 2 .-2(o+c)e k d+d 2 . 
in particular, for b = c and d = c 2 

n (i+i)-ss n = < a - 2 > 

fc ef (2Z+1) fc ef (2Z+1) 

b) Let = a k — v^Az, = a k — a/Az, and Ze£ 

PfS fc = Pf5 fc (a fc , a_ fc ) = (A k} A_ k[ + Xp 2 )(A kl A^ + Xp 2 ) (A.3) 

+ {A k] A ki + A- k] A_ k[ + 2Xp 2 )Xx 2 + AV 

Then 

n PiS k (a k a- k ) TT PfS fc (e" fc Og fc e -" fc 0g_ fc ) 

a^a^fc ' e ->0 11 a l a - k 

fe S^(2Z+l) e<0 fe g^(2Z + l) 

fco>o feo>o 

_ cosh 2 (fV^+VA^a+Apa) -sin 2 (fv/A^+xa)) 

cosh 2 fe k 6 ^ A - 4J 



Proof: We first explain why the definition of the product given in the lemma is the right 
one. 

Consider the quadratic perturbation 

H = Hq + XN = ^ e k ata^a + X? ^2 a k<T a k<r 

kcr kcr 

One finds by explicit computation of the trace 

7(\ R\- Tre-^ H Q + XN '> _ jj l+e-^k+A) 

Tre -f3H —11 1 + e -/3e k 

On the other hand, by explicit summation of the perturbation series or equivalently, by 
performing the fermionic functional integral, one obtains 

Z(X,/3) = J e^^^ k °dpc 

U ik^Z e e U dip k<J dip k<T 

ka £k J ka 



k 

ikp — e k — A 



ssC 1 --) (A - 6) 
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In the perturbation expansion the product in the last equation results from 

n ( i — h-\ = e E feCT log ( 1_ ^) = e ~ E r =i V Y. k <r ^ 

The term for r = 1 requires special attention since ifc _ ek is not absolutely convergent. 
This sum enters the perturbation series as the propagator 

O [Xq — X Q , X — X ) — Tre -f 3 H 

evaluated at xq — x' = 0. Here Tip(xox.)ip + (x' x.') is defined to be ip(xox.)ip + (x' x.') if 
xo > x' and — ^ + {x'qx')i^{xqx.) if xq < x' . By explicit computation one finds that for 

x ^ x'o 

C(x - x' , x-*0 = ££ e ik(x-x') e -e t (x -4) [^(-e k )tf - * ) - M^k^o - x' Q )} 

k 



iW(x — x') — ik (x —x' n ) 1 

ik —e k 



k,fcn 



and for = £q 



C(0,x-x') = J,^e^ x - x ')^(-e k ) = lim ^ J] ^(x-x'j-i*^ 



ek 

<=<0 h,k 



-0v 



Here 6* is a step function being 1 for positive arguments and Op(—v) = is an 

approximate step function. Therefore the /co-sums have to be evaluated according to 

y IT— = y e -F^~ ■= lim e~ lk ^-^- = /%(-e k ) = (3 jf^r 

/ j !K -e k / j ifc -e k e ^ ik () -e k < PV r \+ e pe k 

fc k ^<o 

We emphasize this point because the following computation leads to a wrong result 

E l = \^ sym 1 . = iy [ 1 , 1 1 = \^ -e k _ l-e-^ e k 

ik -e k / j ifco-e k ' 2 / j \_ik -e k -ik -e k j / ; k l+p'^ 2 T+e~^k 

k k k k 

Therefore the product in (A. 5) is found to be 

k V lk ° ek / k V k / 

Eoo 1 Y^ s 2/Tn A A 

__ r=l ~ 2-ffc () (ifc -e k ) r gZ^fe ifc -e k ^feQ ifc 0- e k 



3 -#A_ T-r fcg + (e k +A) 2 „_£ A 



= n S9m fl-T^) e"2 A = n "°T^r ; e~*- 



fe o >0 K o+ e k 



(7/7.25) cosh4(e k+ A) g _ |A = l+e^) (A . 6) 
cosh|e k 1+e pe k v 7 
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which proves (A. 2). 

In the general case one may proceed similarly to obtain (using a k + ci- k = — 2ek) 



jj a k a- k +ba k +ca_ k +d 



lim n ( 1 + »°"*°°»+~-" fc °°-»+'* > ) =: n e 



£<0 

fe a,_ k Z— 'fc ^— 'fcr) afc 

fc 

(a fc a„ fc +ba fc +ca„ fc +d)(a fc a_ fc +ba_ fc +ca fc +rf) ^(b+c) 



T-r y^fc"-— fc I ^fc I fc i Lt.y V ^ fcLt ._ fe | i i it.; _ 

fc >0 Ka-fe) 2 

(afca-fc) 2 +afca- fc (-2e k (b+c)+2d+b 2 +c 2 )+bc(ag+a!. fc )-2e k d(b+e)+d 2 | (6+c) 

Since a 2 . + a 2 fc = — 2/cq + 2e\ = —2a k a- k + 4e k one obtains 

■pre a fc q_ fc +ba fc +ca_ fc +d _ „ (a fc a_ fc ) 2 +pa fc a_ fc +g £(b+ c ) 

fc 6f(2Z+l) afea " fe " ~ fco>0 ~ ( a * a -*) 2 



= n fl + ^-Ul + ^— ) e^ b+c ) 



k >0 

where p = -2e k (6 + c) + 2d + (b + c) 2 - Abe and g = 46CC 2 , - 2e^d(b + c) + d 2 . Application 
of (III. 25) proves part (a) of the lemma. 

Part b) First, if bk = a k — \/Xv, one computes that 

PfS k = (b k b- k + \p 2 ) 2 + {bl + b 2 _ k + 2Xp 2 )X(w 2 + x 2 ) + X 2 (w 2 + x 2 ) 2 (A.7) 
Namely, since 

A k ^A-ki = cikci-k - VXv(a k + a-k) + i^Xw{a k - a_ fc ) + X(v 2 + w 2 ) =>- 



{A k] A_ k[ + Xp 2 )(A kl A_ k} + Xp 2 ) = 



a k a- k - VXv(a k + a-k) + X(v 2 + w 2 + p 2 ) 
+ Xw 2 (a k - a_ fc ) 2 
[b k b- k + X(w 2 + p 2 )} 2 + Xw 2 (b k - b- k f 



and 

A k] A k[ = (a k - VXv) 2 + Xw 2 = b 2 k + Xw 2 

one has 

PfS k = [b k b- k + X(w 2 + p 2 )] 2 + Xw 2 (b k - b- k ) 2 + [b 2 k + b 2 _ k + 2X(w 2 + p 2 )] Xx 2 + AV 
= [b k b_ k + Xp 2 } 2 + (b 2 + b 2 _ k + 2Xp 2 )X(w 2 + x 2 ) + X 2 (w 2 + x 2 ) 2 
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Now let y 2 = w 2 + x 2 and 

F k (a k , a_ fc ) = [b k b. k + Xp 2 } 2 + (b 2 k + b 2 _ k + 2Xp 2 )Xy 2 + X 2 y 4 
Then part (b) follows from the following formula 



Um TT F fc ( e " fc 0afc, e -' £fc 0q^ fc ) = cosh 2 (|V (e k +V\v)*+\p* ) -sin 2 (j ygg) 

11 a fe«-fc cosh 2 f e k ' 1 ' J 

e>0 fe ef (2z+i) 

fe >0 

Proof of (A. 8): One has A k A_ k = a k a- k — \^Xv(a k + a- k ) + Xv 2 and 

A\ + A 2 _ k = a 2 k + a 2 _ k - 2\fXv{a k + a_ fc ) + 2Xv 2 = 2A k A_ k + (a k - a_ fc ) 2 
which gives 

F k = (A k A_ k + Xp 2 + Xy 2 ) 2 + (a k - a_ k ) 2 Xy 2 

= [A k A- k + Xp 2 + Xy 2 + iX^y(a k - a_ fc )] [A k A- k + Xp 2 + Xy 2 - iX^y(a k - a_ fc )] 
= [a k a-k - (VXv - iX^y)a k - (VXv + iX^y)a- k + Xv 2 + Xp 2 + Xy 2 ] x 
[a k a- k - (\f~Xv + iX^y)a k - (\^Xv - iX^y)ci- k + Xv 2 + Xp 2 + Xy 2 ] 

Therefore one obtains 

e F k {a k ,a_ k ) _ ^. a k a_ k -(VX v-iX 2 y)e ik o e a k -(VX v+iX 2 y)e- ik o e a_ k +Xv 2 +Xp 2 +Xy 2 

k >0 < a -* ~°fcj aka - k 

(A.) cosh | v/e 2 +a. + cosh | ^e 2 e _/3^/\ v 
cosh ^e k cosh l[e k 

where 

p = (e k + 2y/\v ) 2 - e\- 4(Aw 2 + Xy 2 ) + 2Xv 2 + 2Xp 2 + 2Xy 2 

= 4v / Afe k + 2Xv 2 + 2Xp 2 - 2Xy 2 
q = 4(Xv 2 + Xy 2 )el + AVXve k (Xv 2 + Xp 2 + Xy 2 ) + (Xv 2 + Xp 2 + Xy 2 ) 2 

= (2VXve k + Xv 2 + Xp 2 + Xy 2 ) 2 + 4Xy 2 el 

such that 

£-q= -4(2VA ve k + Xv 2 2 + Xp 2 )Xy 2 - 4Ay 2 e 2 = -4 [(e k + VXv) 2 + Xp 2 ] Xy 2 
and therefore 

e k + u± = (e k + VXv) 2 + Xp 2 - Xy 2 ± 2iX^y [(e k + VX v) 2 + Xp 2 ] * 
= {[(e k + v / Az;) 2 + Ap 2 ] i ±zA^} 2 
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which results in 



h§ j[(e k +V^) 2 +Ap 2 ]5-HA3yj coshf |[(e k +v^«) 2 +Ap 2 ]i-iAiy| 



ne F k (a k ,a_ k ) _ ^ J J —8\/\v 
— 77^77^ — TT3 7773 e 

fc >0 a fc a -fe cosh£e k cosh£e k 



sinh 2 (4iv^y) + cosh 2 (f v /( e k + V / A«) 2 +Ap 2 ) -fi^\ v 
cosh 2 
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